Critical numbers of attractive Bose-condensed atoms in asymmetric traps 
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The recent Bose-Einstein condensation of ultracold atoms with attractive interactions led us to 
consider the novel possibility to probe the stability of its ground state in arbitrary three-dimensional 
harmonic traps. We performed a quantitative analysis of the critical number of atoms through a full 
numerical solution of the mean field Gross-Pitaevskii equation. Characteristic limits are obtained 
for reductions from three to two and one dimensions, in perfect cylindrical symmetries as well as in 
deformed ones. 



PACS: 03.75.Fi, 32.80.Pj, 42.65.Tg, 02.60.Jh 

The predicted collapsing behavior of the condensed 
systems with attractive two-body atomic interactions , 
first observed in experiments with ^Li was recently 
tested in experiments with *^Rb |3|. In the experiments 
described in Ref. and more recently in by means 
of Feshbach resonance techniques |^ , the two-body inter- 
action was tuned from positive to negative values. Be- 
sides the fact that the experimental results qualitatively 
agree with the theory, and confirm results of previous 
variational treatments , they also show a 

consistent quantitative deviation of about 20% from the 
mean- field predicted critical number of atoms, A^c [0- 
The asymmetry of the trap was shown in Ref. [ p2| to be 
responsible for about 4% of the observed deviation. 

In this respect, it is relevant to obtain precise and 
reliable numerical results for the mean-field calcula- 
tions, in order to probe their consistency and possi- 
ble limitations. The actual experimental atomic traps 
are in general harmonic and non-symmetric. Ex- 
treme asymmetric traps have been recently employed 
in experimental investigations with condensates con- 
strained to quasi-one (ID) or quasi-two dimen- 
sions (2D) iQ, exploring theoretical analysis considered 
by several authors [ p^|l6[p^ ,|7|, p^ , p9[ . A non symmetric 
three-dimensional (3D) trap is reported in Ref. |2^] , with 
the frequencies given by 2ui = \/2w2 — uj^, — 2-11 x 33IIz. 

Considering the general non-symmetric traps that have 
been employed, the accuracy of the comparison between 
experiments and the results of mean-field approxima- 
tion relies in precise calculations using arbitrary three- 
dimensional traps. In case of attractive two-body inter- 
action, the maximum critical number of atoms for a sta- 
ble system is one of the interesting observables to study, 
which is also related to the collapse of the wave-function 
of the system. In these cases, where the two-body scat- 
tering length is negative and the kinetic energy cannot 
be considered to be a small perturbation, the Gross- 
Pitaevskii mean field approximation has been applied, 
given reliable results in explaining the observations in 
the stable (non-collapsing) conditions [^ . 

Before presenting the mean-field equation for an arbi- 



trary 3D case, let us analyze qualitatively the collapse 
phenomenon for asymmetric traps. The interaction en- 
ergy is proportional to the square of the density, varying 
with the negative two-body scattering length. For traps 
with cylindrical (or almost cylindrical) shapes, there are 
two quite different situations: one, pancake-like, with 
the frequencies in the transverse directions being smaller 
than the frequency in the longitudinal direction; the 
other, cigar-like (quasi-lD), with the frequency in the 
longitudinal direction smaller than the frequencies in the 
perpendicular directions. For a true ID system, one does 
not expect the collapse of the system with increasing 
number of atoms However, it happens that a re- 

alistic ID limit is not a true ID system, with the density 
of particles still increasing due to the strong restoring 
forces in the perpendicular directions |^6| . The relevance 
of the quasi-lD trap have been pointed out in Ref. [^ , 
to control the condensate motion. But, as we are going to 
see, the critical number of particles in the quasi-lD limit 
is smaller than the critical number of particles in the 2D 
limit, if we just exchange the longitudinal and perpendic- 
ular frequencies. The physical reason for that behavior is 
the increase of the average density in the cigar-like con- 
figuration relative to the pancake like one for the same 
number of atoms, implying in a strong collapsing force in 
the first case and consequently the cigar-like geometry is 
a more unstable configuration compared to the pancake 
like one. This conclusion is in apparent contradiction 
with the remark made in section IV of Ref. [ p^ , saying 
that, considering the better collapse-avoiding properties, 
"the cigar-shaped trap is the optimal one" . We are going 
to discuss this problem in detail and clarify this issue. 

In the following, we revise the Gross-Pitaevskii (GP) 
formalism, for an atomic system with arbitrary non- 
spherically symmetric harmonic trap. The Bose-Einstein 
condensate, at zero temperature, in the GP mean-field 
approximation is given by 
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m is the mass of the atom, 



/i is the chemical potential, and the wave-function 5" = 
^(r, t) = ^'(r, 0)exp(— i/Ltt/?i) is normalized to the num- 
ber of particles N. The arbitrary geometry of the trap is 
parameterized by three different frequencies wi, ^2 and 
0^3. For convenience, it is appropriate to define the fre- 
quencies according to their magnitude, such that, in the 
present work we assume uJi < uj2 < i^s- 

Here we will be concerned only with systems that have 
attractive two-body interactions [o = — |a|, in Eq. (p" 
In this case, it was first shown numerically, in Ref. 
that the system becomes unstable if a maximum critical 
number of atoms, iVc, is achieved. We present precise 
results for a critical parameter k, directly related to the 
maximum number of atoms, in a general non-symmetric 
configuration of the trap. 

By rewriting Eq.(nl) in dimensionless units: 
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where r = tot, r,; = Iq Xi and (j) = /oy47r|a|^', with 



(rx\(t)\' = 47r- 
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The oscillator length is defined in terms of uj, which is 
taken as the geometrical mean value of the frequencies: 



Hz. With this symmetry, they have obtained k = 0.459 ± 
0.012 (statistical) ±0.054 (systematic). In Ref. [||, 
assuming the cylindrical symmetry wi = 27r x 6.80Hz 
UJ2 = — 21^ X 17.35 Hz, it was obtained k — 0.55, a 
value about 4% lower than kg. 

In our numerical approach, the calculation is per- 
formed by evolving the nonlinear equation (^ through 
imaginary time |24|. The evolution is performed for an 
initial value of the normalization (^ till the wave function 
relaxes to the ground state. The wave function is renor- 
malized after each time step. The process is repeated 
systematically for larger values of the normalization, un- 
til a critical limit is reached. At this critical limit the 
ground state becomes unstable. The time evolution is 
done with a semi-implicit second order finite difference 
algorithm. An alternating scheme is used in the xi and 
X2 direction, with a split step in the 2:3 direction. This 
procedure is done only for Xi > 0, taking advantage of 
the refiection symmetry of the ground state. We consider 
100'^ grid points and time step equal to At = 0.001, veri- 
fying that the algorithm is stable for long time evolution. 
As we increment the normalization, approaching the crit- 
ical limit, the wave function starts to shrink. So, in order 
to maintain the precision, we introduce an automatic re- 
duction of the grid sizes, Axi, Ax2, and Axs, gauged by 
the respective root mean-square-radius in each direction. 

In Fig. 1 we show our main results for the critical con- 
stant k, covering many different geometries. We plot 
k/kg (and k in the rhs y— axis) as a function of A, which 
is defined by 
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For strong non-symmetric cases, particularly when com- 
paring the two extreme cylindrical-shape geometries, 
LUi 102 << 1^3 (pancake-shape) and uji << ^2 ~ W3 
(cigar-shape), it is expected a noticeable difference be- 
tween the corresponding critical number of particles. 

We define a parameter fc, related to the critical number 
of trapped atoms Nc as in Ref. ra: 
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This parameter is a maximum critical limit for stable so- 
lutions of the dimensionless Eq. (||). It will depend only 
on the ratio of the frequencies of the trap. Within the 
precision given in Ref. |Q, kg — 0.5746, where kg is k 
for spherically symmetric traps. In Ref. |Q, the critical 
number was calculated for a nonsymmetrical geometry, 
but in a case that the frequency ratio is not too far from 
the unity (wi/wj, = 0.72), giving a result for the number 
of atoms almost equal to the spherical one. 

In the experiments with ^'^Rb they have considered 
an almost cylindrical "cigar-type" symmetry, with the 
three frequencies given by 17.47 Hz, 17.24 Hz and 6.80 
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FIG. 1. Critical constant k — Nc\a\{mu! /hy^^ , calculated for an art 
> ijJ2 > uJi- The ratio u)zl^\ is shown below each corresponding ci 
dotted-lines correspond to cylindrical (77 = 1) and deformed cylindrical (r 
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We show several curves in which we kept constant the 
ratio ws/wi (sohd-hnes). The values of uj^/loi are indi- 
cated inside the figure, just below the corresponding plot. 
The dashed and dotted-lines correspond to cylindrical 
[rj — 1) and deformed cylindrical (77 — 1.414) symmetries 
for the trap, where 77 is the deformation parameter, as 
it will be explained. In the left-hand-side (Ihs) of these 
plots we have ^3 = riuj2 » uJi (cigar-shape); and, in the 
right-hand-side, >> uj2 — rjuji (pancake-shape), fc can 
be determined for any symmetry that was not shown, by 
interpolating the already given results. The results for 
the complete non-symmetric case are consistent with the 
previously obtained values in cylindrical symmetry [ p2| . 
The maximum value for the critical number k is obtained 
for the spherically symmetric case (wi =012 = W3). 

As already observed, Fig. 1 includes previous calcula- 
tions in the limit of the quasi ID (cigar-shape) and quasi 
2D (pancake-shape) symmetries (see dashed-line) . How- 
ever, for the sake of comparison with previous results 
obtained by several authors, we also present in Fig. 2 
the cylindrical pancake- type (cj^ — loi = 002 « and 
cigar- type [lo± = W3 = 0^2 >> t^i) results. In Fig. 2, for 
cylindrical geometries, we compare our exact results, for 
= Nc\a\{mijj^/hY/'^ as a function of A, with the corre- 
sponding variational ones ofRcfs. |0,|l^. The variational 
results (dashed-lines) are consistently a bit higher than 
the exact ones. For A ^ (wi « ijJ2 = uj^ — 
the exact and variational results for the critical con- 
stant are Nc\a\{mu}i_/liY/'^ =0.676 and 0.776, respec- 
tively. They are consistent with the quasi-lD limits given 
in Refs. ||l^,|l^. When oj^ » uj2 — uji, the variational 
2D limit \/7r/2 of Ref. is comparable with our exact 
result (fcs = 0.931-y/7r/2). In this case, the quasi-2D limit 
coincide with the true 2D limit l2lJl3,E2,E6| . 



like trap, Nc also increases. We need to clarify this mat- 
ter, which is the subject of the next four paragraphs, 
because such result is apparently contradicting a remark 
made in Ref. jl^, saying that "the cigar-shaped trap is 
the optimal one" . 

As shown in Ref. jlj] and also in the present calcula- 
tion, in a deformed cylindrical symmetry, the cigar shape 
(with one of the frequencies smaller than the other two) 
is more favorable to obtain a larger value of k than the 
pancake-shape symmetry (with one of the frequencies 
larger than the other two). See, for example, in Fig. 1, 
the two extreme points of the curve with uj^lji — 100. 
This results from the definition of k, Eq. (ph, in terms 
of the averaged oscillator length Zq. However, the max- 
imum value of k can only be directly related with the 
maximum value of Nc in case that lo is kept fixed. And, 
with /q fixed, Nc is maximized for A = 1, corresponding 
to the spherically symmetric case {k = ks). If we fix 
and ^3/0;!, Nc is maximized for a deformed cylindrical 
symmetry with A < 1, as one can see from Fig. 1. 

Considering exact cylindrical traps, by exchanging the 
frequencies (which, obviously, does not keep constant the 
averaged frequency w), it is valid the following ratio that 
was obtained in Ref. |12| : 
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This result favors the pancake-like symmetry (A — 
ujs/uji — u!3/uj± > 1), to obtain a larger value for 
Nc- Consider, for example, a cylindrical pancake-type 
trap with A = 100, in comparison with a cigar- type 
trap with A = 1/100. We notice that, in this case, 
i?(A = 100) « 1.6, implying that with such pancake-like 
trap (ix'3 = lOOo;^) one can obtain about 60% more par- 
ticles than with the corresponding cigar-like trap (lu^ — 
uj± = lOOwi). Let us consider the recent experiment with 
quasi- ID (cigar- like) trap used in the formation and prop- 
agation of matter wave solitons, with ^Li In this case, 
it was used axial and radial frequencies, respectively, 
equal to 3.2 Hz and ~625 Hz or, uj± — LU3 = 27rx 
625 Hz, uji = 2ttx 3.2 Hz and A = oji/ujs = 0.00513. 
So, as shown in Fig. 2, we are practically in the limit 
A = 0, which gives N^^\a\/ y/h/ (mu;^) w 0.675. Consid- 
ering that the scattering length was tuned to a = — 3ao 
(flg is the Bohr radius), the maximum number of atoms 
in this quasi-lD trap is N^^ w 6400. If we exchange 
the radial and axial frequencies in this experiment, going 
from a cigar-like to a pancake-like trap, uj± = 27rx 3.2 
Hz, UJ3 — 27rx 625 Hz and A = ut^/uj^ = 195.3. In this 
case, Nc'^\a\/ y/h/ {muj-s) « 1.12. So, the critical number 
Nc\a\{muj3/hy^^, for cylindricailntit^ 4^B3^^B3i)rftite isoBjpOT4d6#%hlffiEig0iS]bh 
In the 2D pancake limit, u;± ^e^snd&igr^iufeberfef^'atoffigldiPtllte (piaM- 
Note that A = 0)3/10-1, when 0)2 ^siOQCitei when UJ2 ~ i^a- 

The discussion about the best way to distribute the fre- 
Considering the present analysis, we observe that Nc quencies to obtain the maximum number of atoms was 
in a cigar-like (quasi-lD) trap is smaller than Nc in a first considered in Ref. arriving that the best is to 
pancake-like (quasi-2D) trap. And, as we deform a cigar- do a spherical trap, for achieving maximum density. (In 




100 



FIG. 2. Exact results of fca 
tional approaches (dashed-lines) 
u)± = u)2 = and fca — kX 
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their almost spherical trap they have obtained from ^--^600 
to ~1300 atoms, in an overall agreement with theoretical 
predictions The point is that, to say what is the 

best configuration for the maximum number of atoms we 
must also say what is the constraint in the frequencies. If 
the constraint is to keep fixed the product luiU!2'-^3 than 
the maximum number will be given for the maximum k 
[as shown in Eq.(|^)], that happens in a spherical symme- 
try, in agreement with Ref. But, if two frequencies 
are equal and given a priori, than the best configuration 
to increase N^. is to make the third frequency equal to 
zero, in agreement with Ref. JlGf . However, if only one 
frequency is kept fixed, than the extreme pancake shape 
will contain more atoms, by making the remaining fre- 
quencies go to zero. 

As a final remark about the maximization of TVc, we 
note that, for attractive two-body interactions, the fre- 
quencies in all three directions cannot be initially zero. 
In the experimental process of condensation, the atoms 
must be trapped while going from higher temperature 
T to the critical one, Tc- By making the frequencies too 
small, the Tc will be smaller and the condensation harder 
to achieve [as shown in Ref. [^^, Tc ^ hLu{Ntot)^^^ , where 
Ntot is the total number of particles in the trap] . Once 
the condensate is achieved, we can switch ofi' some of the 
frequencies, and obtain a finite maximum value for Nc, 
as observed in Fig. 2 for T = 0. In the present work, we 
are concerned with the case of zero temperature. Effects 
of temperature on the collapse and on Nc, for spheri- 
cally symmetric case, have been discussed, for example, 
in Ref. [ p8| . Nc is maximized by making uj as small as 
possible. However, a discontinuity exists when taking 
all the three frequencies exactly zero. Without trap the 
collapse will occur for any number of particles. 

We must stress that Fig. 1 also includes novel non- 
symmetric cigar- type geometries ujs = r]uj2 » i^i and el- 
liptical pancake-type geometries uj^ >> uj2 ^ V^i , where 
77 is the excentricity of the ellipsis. In particular, the dot- 
ted line in Fig.l gives the critical values of k for 77 = \/2. 
Thus, one could replot, as in Fig. 2, a series of curves 
whose limiting cases describe non-symmetric cigar- or 
pancake-like symmetries. In the ID and 2D limits, one 
can obtain critical quasi-lD and quasi-2D numbers, for 
each rj. 

One can also obtain an interesting result, for asymmet- 
ric cigar-shape traps, where W3 = r/W2 »Jih-, following 
the approach given in section II of Ref. The as- 

sumption made in [ p^ to reduce the 3D equation to a ID 
equation can be justified in the limit when the forces in 
the transverse directions are the main responsible for the 
trapping potential. So, in the limit loi — s- 0, we are able 
to generalize the solitonic ID equation for a cigar-shaped 
trap, deformed by a given ratio rj = 0^3/^2 between the 
transverse direction frequencies. In this case, rj > 1 can- 
not be arbitrarily larger. Also, by a scaling procedure 
applied to the equation for different values of 77, using 
our previous exact result, we have 
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0.676. (8) 



This generalizes the cigar-shape quasi-lD results of jl^] 
(with r\ such that we still have uj2 >> wi). The ID 
solid line given in Fig. 2 is also applied to deformed cigar- 
shaped symmetries if we replace fcs by fcsT^"^/^ in the 
?/— axis. This result may be relevant for asymmetric 
waveguide propagation as one can deform the cigar-type 
symmetry and control the collapsing condition. From 
Eq. (^) we observe that the maximum critical number 
Nc will increase when deforming the cigar-like symmetry 
by a factor proportional to {^3/102)^^^- 

Now, let us see the effect of deformation in the same 
example of the cigar-shaped geometry used in Ref. [||, 
with ''Li gas. We use the same value of a = — Soq, with 
[^3 — 27rx 625 s""'^, and take UJ2 = ^s/r]- As given in 
Eq. (|), we go from Nc = 6392 (77 = 1) to Nc = 6970 
(77 = V^), or Nc = 7601 (77 = 2). It means, an in- 
crease of ~9% when 77 = -\/2; and ~19% when rj — 2. In 
case, 7] >> 1, the approximation considered in the wave- 
function separation, as given in |T^, is not valid. In 
such case we are reaching the other deformed pancake- 
like symmetry, where lu^ » lu2 ^ uJi. However, in the 
pancake geometry, the effect of deformation in Nc is neg- 
ligible. By comparing the quasi pancake-like geometry 
with the quasi cigar-like geometry, the number Nc in the 
cigar-like geometry is much more sensitive to deforma- 
tions. 
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FIG. 3. The three components, i = 1,2, 3, of the mean-square-radius, 
N\a\/lo, for the case that lus — \/2 uj, L02 ~ oJ, and uii = . 

We have also studied the behavior of the root mean 
square radius for the case 0^3 = '\/2w2 = 2<x'i. We verified 
that, as the system approaches the critical point (or col- 
lapse), the wave function tends to be more "spherical", 
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confirming earlier conclusion made with a Gaussian vari- 
ational approximation |3C[| . In Fig. 3, we show the corre- 
sponding results, for the three components of the mean- 
square-radius, (rf) = Iq{x1) [i ~ 1, 2, 3). As shown, when 
iV = 0, we have {r\) / {r^) = {rl)/{rl) = \/2; and, when 
N fa Nc such ratio is drastically reduced. 

In summary, we have calculated systematically the 
critical number of particles, in systems that have neg- 
ative two-body interactions, for traps with arbitrary ge- 
ometries. The maximum critical number of particles Nc 
can be derived from the given value of the parameter 
k, given in Eq. (||), once one has the scattering length, 
the mass of the atomic system, and the frequencies of 
the trap. The results are shown in Figs. 1 and 2. The 
value of k, for any symmetry non explicitly given, can 
be easily derived from Fig. 1, by interpolation. It is 
also pointed out that the results shown in Fig. 2, in the 
ID cigar-like case, can be extended for slightly deformed 
cylindrical symmetries, by replacing the axis label k^ 
by ksT]-^/*. 

Our main results in the present work are: (i) The 
maximum number of particles, Nc, for arbitrary 3D trap 
geometries is given through the results shown in Fig.l. 
(ii) The optimal trap configuration, to avoid the collapse 
with maximum Nc, is found to be strongly dependent 
on the constraints of the frequencies of the trap. If we 
initially fix one of the frequencies, the best configuration 
of the trap is pancake-like, with the other two frequen- 
cies going to zero. Analogously, if we initially fix two 
equal frequencies, the best configuration of the trap is 
cigar-like, with the third frequency close to zero. If we 
initially fix two different frequencies and try to vary the 
third frequency between the fixed ones, the best config- 
uration is again pancake-like. We show that Nc is much 
more sensitive to deformations of the trap in a cigar- 
like geometry than in a pancake-like geometry. Finally, 
for small deformations rj of the cigar-like traps, where 
rj > 1, releasing the longitudinal direction, the solitonic 
solutions, obtained in Ref. [|6| will be rescaled by the 
deformation. Nc will be rescaled by a factor Ty^/**, gener- 
alizing the findings of Ref. jl^ . 
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